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The calibration of the in-flight antenna beam shape and possible beam degradation is one of 
the most crucial tasks for the upcoming Planck mission. We examine several effects which could 
significantly influence the in-flight main beam calibration using planet transit: the problems of the 
variability of the Jupiter’s flux, the antenna temperature and passing of the planets through the 
main beam. We estimate these effects on the antenna beam shape calibration and calculate the 
limits on the main beam and far sidelobe measurements, using observations of Jupiter and Saturn. 
We also discuss possible effects of degradation of the mirror surfaces and specify corresponding 
parameters which can help us to determine these effects. 


1. INTRODUCTION 

The ESA Planck Surveyor 0 is designed to im¬ 
age the whole sky of the Cosmic Microwave Back¬ 
ground (CMB) radiation with unprecedented sensi¬ 
tivity (A T/T ~ 2 x 10 -6 ) and angular resolution 
(down to 5') at 9 frequencies: 30, 44, 70, 100, 143, 
217, 353, 545, 857 GHz at Low Frequency Instrument 
(LFI) p] and High Frequency Instrument (HFI) [2]. 

To achieve these high sensitivity and resolution, 
it is necessary to carefully account for all potential 
systematic features in the data [3]. 

One of the systematic effects is related to the 
in-flight antenna beam shape and its reconstruction. 
Apart from the need to acquire the radiation pattern 
of the antenna beam before the flight, the calibration 
of the in-flight antenna beam shape is one of the key 
components for achieving the scientific goals of the 
Planck mission. This problem has been considered by 
different Planck groups H 0 0 H|, and recently by 
01 - 

The accuracy of the CMB anisotropy Ci estima¬ 
tion will be affected, among other experiment param¬ 
eters, by our ignorance of the in-flight antenna beam 
shape of the main beam and far sidelobes, and pos¬ 
sible degradation of the mirror surface shapes. When 
scanning the antenna beam of the Planck mission is to 
move across the sky, implying that the antenna beam 
response is a function of time. After pixelization of 


1 http://astro.estec.esa.nl/SA-general/Projects/ 
Planck/ 


the time-ordered data the position of each pixel on 
the pixelized CMB map is directly related to some 
data points in the time stream. It is therefore nec¬ 
essary to obtain the information about the in-flight 
beam shape, its inclination and the location of the 
beam center relative to each pixel, in order to im¬ 
prove the model of the in-flight main beam shape as 
well as in the far sidelobe region. 

To tackle this issue, Burigana et al. [6] have sug¬ 
gested a method which uses planet transit to recon¬ 
struct the in-flight beam shape. These planet cross¬ 
ings can help the in-flight beam recovery down to 
—25-732.5dB at 30GHz. They also showed that the 
main beam pattern can be described by the bivariate 
Gaussian approximation. Recently Chiang et al. [7] 
have proposed another method for the beam shape 
estimation based on the interplay of amplitudes and 
phases of the sky signal and instrumental noise. This 
method is useful in extraction of the antenna main 
beam shape down to —7-710 dB, and does not need a 
strong radio source shape calibration. These methods 
have laid a base for the determination of both the in¬ 
flight antenna beam shape and its variations during 
observations. 

The aim of this paper is to re-examine in details 
the proposed method of the in-flight antenna shape 
reconstruction by planet crossing 0 in order to esti¬ 
mate possible beam degradation effects. The method, 
based on Jupiter and Saturn observations, has some 
subtleties needed to be addressed such as the tempo¬ 
ral and frequency variations of the Jupiter’s flux and 
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Saturn 



Figure 1: Positions of Jupiter and, Saturn at the be¬ 
ginning of the Planck mission (November, 2008). The 
radius of each concentric circle is differed by 2 AU. 

passing of the planets through the main beam. 

This paper is arranged as follows. In Section 2 we 
describe the Planck scan strategy in relation to planet 
observations, and the beam definitions. In Section 3 
we discuss the in-flight beam calibration using planet 
transit. We concentrate in Section 4 on the effects in 
beam calibration due to fluctuation of the Jupiter’s 
flux density. In Section 5 we discuss variations arising 
from the distance to the planets and from the scan 
strategy. Conclusion is in Section 6. 

2. SCAN STRATEGY AND THE 
BEAM PROBLEM 

The proposed scan strategy [4] for the Planck mis¬ 
sion is a whole-sky scanning with the satellite orbit¬ 
ing around the L2 Lagrangian point of the Earth-Sun 
system. The satellite spin axis will be pointed in the 
anti-Sun direction and will have a tilt to the ecliptic 
plane of 10°. The telescope optical angle is inclined 
at 85° to the spin axis. The telescope will scan the 
same circle 60 times around the spin axis at 1 r.p.m. 
Each hour the spin axis will be shifted along the eclip¬ 
tic plane with a step of 215 during the entire mission 
(about 15 months). 

2.1. Observations of the Planets 

According to the present schedule, the launch of the 
satellite is planned in June 30, 2008 and the flight to 
the L2 point will take approximately six months. To 


construct the L2, Earth, Jupiter and Saturn orbiting 
around the Sun, we use the planets ephemerides cal¬ 
culation procedure of the XEphem.3.7.2 software by 
[TO] (see Fig. [T|. 

As mentioned above, the scan angle is assigned 
to be 85°. Assuming the beginning of the mission 
on the 31th of October 2008, in 15 months Jupiter 
will cross the main beam direction 3 times on the 
dates of 22.05.2009, 03.11.2009 and 27.01.2010 (and 
on 09.12.2010) with an accuracy of ±1.5 day, mainly 
due to tilt and starting conditions, whereas Saturn 
will cross the main beam 3 times on the dates of 
17.12.2008, 31.05.2008 and 30.12.2009 (and in 2010, 
on 11.06.2010) (Fig.®. 

As shown in Fig.0 the antenna beam shape cali¬ 
brations using Jupiter and Saturn can only be realized 
6 (3 + 3) times during the mission with a few differ¬ 
ent time intervals between each planet crossing. These 
“windows” provide us with the possibility of estimat¬ 
ing degradation of the antenna main beam shape from 
the long time frequencies (up to 15 months) to the 
short frequencies (about 10 days)H For continuous 
calibration of the main beam area, therefore, it is 
necessary to use the method by [7] to reconstruct its 
ellipticity and orientation in the data analysis. One 
of the interesting dates in the planet observation is 
03.11.2009, when Jupiter and Mars could be observed 
during one day. Unfortunately, the unpredictability 
of the Mars spectrum in the millimeter wavelength 
range due to atmosphere storms m does not allow 
us to calibrate the far sidelobes. However, the mission 
can provide interesting results for the study of Mars 

m- 

2.2. Inclination of the Planet Orbits to the 
Ecliptic Plane 

We make our estimation of the planet transit for the 
ecliptic plane projection. However, the planets have 
ecliptic latitude different from zero. The maximum 
inclinations are « 1° for Jupiter and « 2° for Saturn 
(the detailed calculation for each planet crossing is in 
Table Q]). Here we note that this problem is not essen¬ 
tial for our case but is similar to the tilt projection 
problem discussed below. 

According to the scan strategy proposals the spin 
axis can have a tilt to the ecliptic plane of about 
10°. This means that the usual ecliptic projection of 
the scan angle will be narrower and we shall observe 


2 The interesting point of the track crossing between Jupiter 
and Saturn near the scan angle on 26.05.2009 allows us to dis¬ 
cuss the possibility of simultaneous observations of two planets 
in one scan. However, to make such observations, we need to 
have a scan angle of about 80°. Otherwise, these transits can 
be used to test the far sidelobes of the beam. 
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Figure 2: Transit of Saturn, Jupiter and Mars through a scan angle of 85°. The dashed, solid and dotted curves 
represent the Saturn, Jupiter and Mars paths, respectively. The solid horizontal line demonstrates a scan angle 
of 85°. Dates near the horizontal line correspond to the Jupiter’s and Saturn’s crossings of the 85°scan angle. 


Table 1: Position of the planets on the dates of crossing by Planck on the ecliptic plane and antenna tempera¬ 
tures (T a = flpianet/^beam • Tb) in transit via scan angle of 85° 


Planet 

Date 

85° 

Size 

(arcsec) 

30 

33' 

T a (mK’ 
100 i 
10' 

1 

>45 GHz 
5' 

Ecliptic 

latitude 

Ecliptic 

longitude 

Saturn 

17.12.2008 

17.9 

7.534 

86.98 

291.2 

+1:58.0 

171:27.2 

Jupiter 

22.05.2009 

40.4 

43.86 

556.3 

1672 

-0:45.6 

325:59.8 

Saturn 

31.05.2009 

17.9 

7.534 

86.98 

291.2 

+2:07.8 

164:57.5 

Jupiter 

03.11.2009 

40.9 

44.40 

563.2 

1693 

-1:01.2 

317:46.9 

Saturn 

30.12.2009 

17.7 

7.450 

86.01 

287.9 

+2:16.9 

184:19.6 

Saturn 

17.06.2010 

17.5 

7.366 

85.04 

284.7 

+2:21.9 

177:57.1 

Jupiter 

27.06.2010 

40.4 

43.86 

556.3 

1672 

-1:16:03 

2:10:53 


Jupiter on other days with another antenna temper¬ 
ature. The scan projection ip of the scan angle a for 
the tilt angle </> can be calculated as 

tan 0 = tan <p cos a, (1) 

which gives us ip = 84.92°, for a = 85° and </> = 
10°. This difference is within 5', indicating that the 
projections of the planet trajectories and the tilt of 
the axis will not interfere with our estimates. 

2.3. Beam Descriptions and its Variations Due 
to Mirror Degradation 

The Planck “antenna beam” is usually referred to 
as a physical model of the antenna response and its 
ground-based verification before the flight, and the in¬ 
flight antenna beam as the beam reconstructed during 




Figure 3: Beam degradation parameters. The left, 
panel illustrates the orientation angle of the beam, 
and the right panel shows cr_ and cr+, the semi-minor 
and semi-major axes of the elliptical main beam. 
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North Ecliptic Pole 



Figure 4: The Planck scan strategy. The square A is 
located at the parallel scan area, where the planets can 
be found. The square B has crossing(s) of the circular 
scans. 


to conduct a detailed investigation of the correspond¬ 
ing degradation of the antenna beam shapes, using 
Jupiter’s and Saturn’s transits as suggested by [6) 3 I . 
In general, beam degradation can be described by 3 
parameters as illustrated in Fig. [3] First of all, the 
mean beam width is a function of time, 

+ct+) = a 2 = cr 2 {t), (6) 

where cr_ and <j _|_ are the minor and major axes of the 
elliptical main beam. The orientation angle a between 
the scan direction and the major axis of the elliptical 
main beam can also be a function of time, i.e., 

a = aft). (7) 

and so can the ellipticity ratio of the beam 

^ = p = p{t). (8) 

O- 

We can use these 3 parameters as the indicators of 
the degradation level of the in-flight antenna beam. 


the flight, which is crucial for possible beam degrada¬ 
tion estimation. The in-flight antenna beam plays a 
significant role in the Ct estimation as well as in the 
extra-galactic point source extraction |12| . 

Physics optics calculations have shown that the 
main beam is roughly elliptical 00 , so we can ap¬ 
proximate the antenna pattern as a bivariate Gaus¬ 
sian beam. The geometrical property of the beam in 
the time domain can be described as follows. We de¬ 
note by xo and yo the position of the beam in a co¬ 
ordinate system fixed to the detector with x in scan 
direction and y perpendicular to x and the beam axis. 
Then the beam shape can be written as 


B t (x 


x t ) = exp 


-i(RU) T D^(RU) 


with 


( 2 ) 


U = 


X - Xo \ 
2/-2/o ) 


( 3 ) 


where R is the rotation matrix which describes the 
inclination of the elliptical beam, 


R = 


cos a 
— sin a 


sin a 
cos a 


( 4 ) 


with a being the orientation angle between the x axis 
and the major axis of the ellipse. The D matrix de¬ 
notes the beam width along the ellipse major axis, 
which can be expressed as 


D = 



( 5 ) 


If, during the routine operation, the mirror sur¬ 
faces are slightly perturbed (deformed), it is necessary 


3. PLANET TRANSITS AND THE 
PIXEL DOMAIN 


Following [6], we can specify the in-flight Planck an¬ 
tenna beam shape model by using Jupiter as a “stan¬ 
dard candle” for calibration. For the Planck (both 
LFI and HFI) frequency ranges we can model the 
Jupiter’s contribution to the resulting A T(r) sky 
temperature in some direction r as 


ATj(r) 


Sj(V,t) 

( hc V 

2 sinh(r 7 / 2 ) 

2k 

^ AjTcmb ) 

ri 2 




(9) 


where r and r j are the unit vectors in the correspond¬ 
ing direction on the sky and Jupiter’s location in a 
given coordinate system, respectively, 6(rj,r) is the 
Dirac delta function, Sj is the Jupiter’s flux, Tcmb, h, 
k and c are the CMB temperature, Planck constant, 
Boltzmann constant, and speed of light, respectively, 
and r] = hv/kTcMB- 

Each observed time-ordered subscan m\ that in¬ 
cludes Jupiter’s image is related to ATj(r) from 
Eq. ([21) through a convolution with the antenna beam 
function B{r,r ) (see [7]) 


yi = df t + n\ 


( 10 ) 


where 


d\ = ATj[r(t)]®B[r(t l ),r(t)\, (11) 

where 0 denotes convolution, and n\ now is the CMB 
signal plus all the foregrounds and the instrumental 


3 As we will see in the next section, the Jupiter’s flux fluc¬ 
tuation could mimic the in-flight beam degradation effect as 
well. 
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noise. The index i marks the i-th subscan with the 
same orientation of the spin axis of the satellite. 

Using all *6 [1, N] subscans, where N is the total 
number of the subscans, we can define the circular 
scan as some linear transformation of the d\: 

dt= Ad\ (12) 

where d* = {dj} is the data vector, and A is the ma¬ 
trix of the transformation. Eq. ED gives us the rela¬ 
tion of a single circular scan for a fixed orientation of 
the spin axis. For simple summation of the subscans 
we obtain d t = (l/A) Below we shall use the 

circular scan as a basic element for the map-making 
algorithm taking into account that the variance of the 
instrumental noise for such a scan is expected to be 
~ iV -1 / 2 times smaller than for each subscan, if the 
instrumental noise is pure white noise. For all circu¬ 
lar scans we can define the vector of the time-ordered 
data y = M s -I- n, where M is the corresponding map¬ 
making matrix, s denotes the pixelized map and n is 
the noise vector [14] , It is worth noting that for the 
in-flight antenna beam shape reconstruction by us¬ 
ing Jupiter’s and Saturn’s images, we do not need to 
construct whole-sky maps, because the — 40 dB limit 
of the expected Planck antenna beam shape corre¬ 
sponds to an angular scale 0f B ~ 5 degrees at 30 GHz 
LFI, 0 and for that purpose we can use the flat-sky 
approximation centered around t he Jupiter’s image 
(Fig. 0]), and apply the method by [7J. Furthermore, 
this assumption allows us to use a linear map-making 
algorithm (see [14)). which is similar, for example, to 
the COBE pixelization schemed 

The signal in each pixel of the map s is then [14] 

s = Wy, (13) 


where W is the corresponding matrix, which depends 
on the scan strategy of the Planck experiment. For 
example, for the simplest COBE pixelization we can 
use W = [M t N _1 M] 1 M T N _1 , where IM = (nn T ) is 
the noise covariance matrix. 0 

Let us go back to the single circular scan. As seen 
from Eq. © and Eq. ED- for the simple average of 
the subscans (A —> 1/A r JA) the Jupiter’s image after 
beam convolution in a circular scan should be 


dt = 


1 


2 kN 


2 smh(j) / 2)hc 


kTcMPjf 


EE 4 5 6 j(Mi)iM0O,(14) 


4 For LFI+HFI frequency range this scale corresponds to the 
minimum. 

5 Note that for estimations of the in-flight beam distortions 
caused by 1// noise, foreground contaminations and so on, [6] 
have to use the whole-sky map for the 1// noise removal. 

6 Note that we can use any modification of pixelization without 
loss of information, e.g. GLESP Q]3[ng or HEALPix ITT 



Figure 6: The power spectrum of the flux density from 
Jupiter at the 13 cm wavelength, which is produced by 
linear interpolation in the parts of the intervals in 
Fig. [H where data are not available. The fitted curve 
is described by a power law P(k) oc k~ 2A . 


where 0, = r (tf)— r(tj,i), is the Jupiter’s location 
in the i-th subscan and B t denotes the beam shape 
in each subscan. As one can see from Eq. in and 
Eq. ®, in the pixelized map the pixels containing 
Jupiter’s image are related to the B t and Sj ( v. tf) and 
can be denoted as follows, 

s oc WS'j(^,t i )B t (0 l ). (15) 

Therefore, for the definition of the antenna beam 
shape in the pixel domain, we can specify some pos¬ 
sible sources of uncertainties from Eq. OSD- 

Before focusing on the two categories of variabili¬ 
ties in beam calibration in the next two sections, we 
would like to briefly mention the uncertainty which 
is related to the location of Jupiter, 


r(tj,i) = + Ar(tj,j) = r(t + S r ), (16) 

where r(tj ) i) indicates the average location and 
Ar(tj,j) corresponds to the fluctuation of the 
Jupiter’s location. Generally speaking, this source is 
related to the pointing accuracy of the Planck experi¬ 
ment. It is natural to assume that (Ar(ij *)} =0, but 
<|Ar(f J , i )| 2 >^0. 

4. VARIABILITY OF THE 
JUPITER’S FLUX ON THE BEAM 
CALIBRATION 

The fluctuations of the Jupiter’s flux can be crucial 
for the in-flight antenna beam shape reconstruction 
scheme. The temporal variations in the Jupiter’s flux 
can be expressed as the constant flux Sj(v), and a 
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Figure 5: Variability of the flux density of synchrotron radiation at f.O) AU from Jupiter at 13 cm wavelength 
(2.3 GHz) (re-produced from fl9 fl. 


fluctuating part ASj(v,t), 

Sj(v,t) = SjO) + ASj(v,t) = «Sj(V)[l + 5 s (t)}. (17) 

Returning to Eq. m and Eq. m in order to 
define the beam shape properties in the pixel domain, 
we will assume (Ar(tj,i)) = 0 (so that in Eq. (fTHl) 
S r = 0). The pixelized beam can be obtained from the 
subscans including Jupiter’s image as follows m- 

B P h) = E E E W [ X + S s(?> (18) 

tep i t ' 

where 7 is the angle between the position of pixel 
which corresponds to the Jupiter’s location in a map 
and position of each different pixel. Possible variation 
of the Jupiter’s flux produces an additional source of 
peculiarities in the pixelized beam shape definition 
proportional to 1 + 6s(t). 

4.1. Characteristic Time Scales 

There are 3 characteristic time scales related to the 
Planck scan strategy. For each subscan of 1 r.p.rn. 
the time scale is T su b — 1 minute. The next time 
scale is that for a circular scan T c i r = 60 minutes, 
which is the time interval for data accumulation in 
one circular scan with a fixed orientation of the spin 
axis. In terms of order of magnitude, T c j r scale can 
be used for estimation of the characteristic time scale 
for the signal variation in one pixel around the main 
beam area. Another time scale is related to the scale 
of the far sidelobes Tps — 0fs/1° days, where 0fs is 
the angular measure subtended by the far sidelobes 
of the beam. For example, for the LFI 30 GHz chan¬ 
nel, the threshold of —30 dB subtends the angular 
scale 0fs — 1-5 degrees [ 6 j and thus Tps — 1.5 days. 
The high-frequency fluctuations of the Jupiter’s flux 
which corresponds to the time scales T su b, T c , r and 
Tps are thus very important for the in-flight antenna 
beam shape reconstruction and may require more de¬ 
tailed investigations, for example by ground-based 


telescopes. These time scales also indicate that all 
irregularities of the Jupiter’s flux for T > Tps cor¬ 
respond to long-term variations and could mimic a 
beam shape degradation effect. 

Unfortunately, we do not know exactly the prop¬ 
erties of temporal variations in the Jupiter’s flux in 
the Planck frequency range 30-857 GHz. According 
to the information available in the literature, the 
Jupiter’s flux variability at a frequency of 2.3 GHz, 
nearest to LFI frequency range, is related to syn¬ 
chrotron emission from the Jupiter’s magnetic belts. 
Fig. H] shows the variation of the flux density of the 
synchrotron radiation from Jupiter at 2.3 GHz, repro¬ 
duced from m- The interval of measurement is one 
day, so any fluctuation shorter than 1 day is yet to be 
measured. In Fig. [ 6 ] we show the power spectrum of 
Fig. El which is produced from linear interpolation in 
the part of the intervals where data are not available. 

The temporal variation that can significantly dis¬ 
tort the beam reconstruction for the 30 GHz channel 
is between 2.4 hours and 1.5 days, the Planck crossing 
time of Jupiter. The variation period of the Jupiter’s 
flux shorter than 2.4 hours may be smeared out af¬ 
ter pixelization in this channel because of the scan 
strategy and beam shape properties. For a variation 
period longer than 1.5 days, the distortion is much 
less. 

In Fig. [7] we show simulated beam reconstructions 
in the 30 GHz LFI channel with possible Jupiter’s flux 
variations. According to Eq. (ED: if 5s(t) is a random 
process, in the Fourier domain it can be character¬ 
ized by a power spectrum. We will assume that for 
Jupiter atmospheric emission the power fluctuations 
5s (t) should have a form P(ui) ~ aw - " + b , where 
a and b are constant. For illustration we will use the 
same values of a and b parameters as for the syn¬ 
chrotron emission discussed above (Fig. ED- The flat¬ 
ness of the power spectrum at large Fourier modes in 
Fig. El allows us to assume the same amplitude of the 
variation for different periods. Panel (a) is the recon- 
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Figure 7: Simulation of beam reconstruction with vari¬ 
ation of Jupiter’s flux. We simulate 30 GHz LFI 
channel (FWHM= 33' with an ellipticity ratio =1.3) 
with (j cmb and, pixel noise cr p i X equal to 3 x 10“ 5 and 
8 x 10^ 6 , respectively. Panel (a) is without Jupiter’s 
flux variation, (b) with Jupiter’s flux varying with a 
period of ).8 hours, (c) 10 hours (the rotation period 
of Jupiter) and (d) 1.5 days (the Planck crossing time 
of Jupiter). The contour lines are annotated in dB. 

structed image without variation in the Jupiter’s flux, 
(b), (c) and (d) are with flux fluctuation of the vari¬ 
ation period equal 4.8 hours, 10 hours, and 1.5 days, 
respectively. The period of 10 hours corresponds to 
that of Jupiter’s rotation. The amplitude of the vari¬ 
ations is assumed to be 20%. We can clearly see that 
the orientation of the main beam changes because of 
the fluctuation of the Jupiter’s flux. 

4.2. Millimeter Spectra of the Planets and 
Their Variations at Different Frequencies 

To estimate the effects of the flux density variation 
of the planets, we have to look at their radio spectra. 
The total Jupiter’s radio spectrum consists of the fol¬ 
lowing two components [21] : a low radio-frequency 
part, which is related with the synchrotron emis¬ 
sion from energetic electrons spiraling in the Jupiter’s 
magnetic field, and a high radio-frequency part, which 
corresponds to the thermal atmospheric emission. 

The synchrotron emission dominates in the fre¬ 
quency range v < 10 GHz while at v > 30 GHz the 
Jupiter’s radio flux is determined by the atmospheric 




Figure 8: The millimeter and sub-millimeter spec¬ 
trum of Jupiter (top) and Saturn (bottom) (repro¬ 
duced from m- The filled diamonds denote the 
fluxes calculated in this work according to the model at 
the 9 observing frequencies in the Planck experiment. 
Note that the first strong dips in the spectra coincide 
around Planck 5)5 GHz observing channel. The figure 
of spectrum was prepared using JiT ]/ (filled circles), 
[If (diamonds), f (squares). Wfj (crosses). 


emission. First of all, let us describe the contribu¬ 
tion of the synchrotron emission and its variability 
to the Ss(tj) parameter for the 30GHz LFI channel. 
Recently simultaneous observations with the Cassini 
spacecraft, the Gallileo spacecraft and the VLA in the 
centimeter wavelength range have been made [18i 119) . 
As is shown by [T5], the Jupiter’s magnetosphere 
is strongly affected by the solar wind. When inter¬ 
planetary shocks propagate from the Sun and reach 
Jupiter, they compress and re-configure the magne¬ 
tosphere, producing a strong magnetic field and elec¬ 
tron acceleration. The ^10% variation with a vari- 
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ability of about 0.5 Jy per month of the 13 cm flux 
from Jupiter observed by PI, and shown in Fig. 0 
may most likely be due to such an effect. This value 
can be used as an upper limit in our estimations. 
For the Planck antenna beam shape reconstruction 
we must, as shown earlier, know the variability of 
the Jupiter’s flux on a time scale of T c j r ~ 1 4- 2 
days, to obtain the lower limit of the variation for the 
Planck frequency range 30-857 GHz. Obviously, the 
synchrotron emission is important, in principle, only 
for 30 GHz channel and it determines the lower limit 
of the Jupiter’s flux variation, if the atmosphere emis¬ 
sion does not produce any fluctuations of the flux in 
the 30 GHz band. We plot in Fig. 0 the power spec¬ 
trum of the 13 cm synchrotron flux variability. From 
this spectrum we find the limit A sync h ~10% per day. 

If we argued that the same amplitude of the syn¬ 
chrotron emission also occurred at the 1 cm wave¬ 
length and used this value in our simulations, it 
would give the lower limit of the flux variation due 
to the synchrotron emission at 30 GHz of 5 sync h ~ 
Cr 6 sy nch/T 6 )A synch ~ 10 -3 , where T b is the bright¬ 
ness temperature corresponding to the total planet 
flux and T'/ ynch that corresponding to the synchrotron 
emission. Thus, we can conclude that the variation of 
the synchrotron emission at v = 30 GHz is not im¬ 
portant for the antenna beam shape reconstruction 
for the whole range of interest (> —60 dB). However, 
it is necessary to obtain additional observational data 
on the intrinsic atmospheric emission. 0 

Detailed studies by m of the millimeter and sub¬ 
millimeter spectra of Jupiter and Saturn have shown 
(see Fig. 0) that there are peculiarities in the spec¬ 
tra in this wavelength range. The two model spectra 
shown in Fig. 0 are from [23], using different physi¬ 
cal parameters such as the size of NH 3 clumps and 
the particle and gas scale heights ratio. Tempera¬ 
tures of Jupiter and Saturn estimated at the corre¬ 
sponding Planck observing frequencies are shown by 
filled diamonds in Fig.0 The first strong dip on both 
spectra almost coincides with the observing frequency 
545 GHz near the 570 GHz NH 3 and PH 3 resonances. 
Estimated brightness temperatures with an accuracy 
of about 10% are given in Table 0 (see also Table 0). 

Unfortunately, we do not have information about 
variability of the Jupiter’s and Saturn’s fluxes in the 
range 30 to 857 GHz, which determines the accuracy 
of the beam shape reconstruction. Some naive expec¬ 
tation of a possible variability in the frequency range 
of interest for the Planck mission could be related 


' We would like to argue that it is natural to expect a non-zero 
fluctuation from the atmospheric emission. Partly our assump¬ 
tion is based on the millimeter and sub-millimeter spectrum of 
the atmospheric emission measurements at 30—857 GHz. 


Table 2: Brightness temperature T b of Jupiter and 
Saturn at Planck observing frequencies with an ap¬ 
proximately 10% accuracy 


(GHz) 

A 

(mm) 

beam 

(arcmin) 

T W er 

Si^urn 

30 

10.00 

33.0 

152 

133 

44 

6.82 

24.0 

158 

135 

70 

4.29 

14.0 

167 

138 

100 

3.00 

10/9.2 

173 

141 

143 

2.10 

7.1 

176 

144 

217 

1.38 

5.0 

179 

146 

353 

0.85 

5.0 

178 

141 

545 

0.55 

5.0 

133 

118 

857 

0.35 

5.0 

145 

114 


with the observed 20% deviations of the Jupiter’s and 
Saturn’s temperature from the pure black body law 
T(i/)=const. For Jupiter, this 20% deviation allows 
us to expect that some process leading to such kind 
of variations can be variable in time at the same level 
and have the characteristic time scale close to the 
period of Jupiter’s rotation (i.e. ~ 10 hours). This 
problem needs an additional and more detailed in¬ 
vestigation by using large ground-based radio tele¬ 
scopes in order to measure possible variation of the 
Jupiter’s and Saturn’s fluxes in the Planck frequency 
range. Recent explorations of Jupiter with the ra¬ 
dio telescope RATAN-600 have shown stability of the 
Jupiter’s flux at the level of 0.1% at the frequency of 
30 GHz for 34 days of observations [26] ■ Such type ob¬ 
servations will be necessary during “in-flight” beam 
calibration of the Planck mission. 


4.3. Expected Polarization of the Flux 


One of the main goals of the Planck mission is the 
CMB polarization measurements. First of all, we 
would like to point out that a low limit of polariza¬ 
tion of the Jupiter’s flux at 30 GHz exists, which is 
related to synchrotron emission. 

According to m, the polarization level n of syn¬ 
chrotron radiation is related to the spectral index [3 

(-Uynch OC I'd') as 


3/? + 
3/3+1’ 


(19) 


which gives ~ 104-75% for different values of n. 
The total polarized flux of Jupiter and Saturn is 
~ Tpianetn. This fact creates a pre-condition to use 
these planets for polarized antenna beam shape cali¬ 
bration. 

Using the value of j3 = —1.26 for v > 13.6 GHz 
for the spectral index of the Jupiter’s synchrotron flux 
pT9] . one can obtain nj up i t er=28%, indicating that the 
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polarized part of the total flux can reach around 0.3% 
in the 1cm wavelength range (see Table 2). 


5. EFFECTS FROM THE STRATEGY 
OF THE OBSERVATIONS 

5.1. Variations of Planet Antenna Tempera¬ 
tures Versus Distance to the Planets 


From the dates of crossing the scan angle by the plan¬ 
ets, one can calculate the corresponding distance and 
hence the angular sizes of the planets. When the ob¬ 
ject size is sufficiently smaller than the solid angle of 
the beam, the antenna temperature T a of the planet 
is given by 


T a = T b 


n 


planet 

^beam 


( 20 ) 


where T b is the brightness temperature of the planet 
(Table [2]) , fl p i an et is the solid angle of the planet in 
observation (angular size in steradian), and flbeam is 
the solid angle of the beam calculated with a simple 
approximation by a Gaussian shape 

B = exp (y^)' <21) 


where a = (al + er+)/2 = d b /2.355 and 0 b is the 
FWHM of the main beam. The results of calculation 
of the antenna temperatures for the three frequencies 
30, 100 and 545 GHz with the corresponding FWHM 
sizes of 33, 10 and 5 arcmin are given in Tabic [T] 


5.2. Scan Strategy and Peculiarities of the In¬ 
flight HFI Beam Reconstruction 

In this subsection we would like to focus on the dis¬ 
cussed scan strategy of the Planck mission and its 
influence on the in-flight antenna beam shape re¬ 
construction using Jupiter’s and Saturn’s transits. 
According to the Planck mission requirement, the 
FWHM for the 10 LFI + HFI channels is shown in 
Table [5J Let us concentrate on the beam shape prop¬ 
erties above — 30 dB for all LFI + HFI channels. 

According to the scan strategy (see Fig.4), the ori¬ 
entation of the telescope spin axis during one hour 
(i.e. 1 r.p.m. of spin for 60 sub-scans) of observa¬ 
tions should be stable: the orientation of the 60th 
sub-scan is parallel to that of the 1st sub-scan at the 
moment to when a given circular scan starts to be 
measured. At the end of the 60th sub-scan, the spin 
axis (and the optical axis) should change its orienta¬ 
tion by 2.'5 in the ecliptic plane (re-pointing). Thus, 
during Jupiter’s and Saturn’s transits the highest res¬ 
olution scale from which the images of Jupiter and 
Saturn can be recovered in the pixelized map is 2.'5 


on one side (as the sizes of both the planets are less 
than 1'). 

Due to the pixelization scheme, however, we 
will face the following two situations concerning the 
higher frequency channels: good and bad cases of the 
planet transit. The good case is where the planets are 
caught by the beam peak just after re-pointing and 
have a maximum signal on a circular scan. The prob¬ 
ability of such a case is low. The bad case is where 
the planets bypass the beam maximum so that the 
point of maximum of the planet flux is missed by 
both of the neighboring circular scans. In terms of 
the map-making algorithm it means that the point of 
the maximum of the Jupiter’s flux is formally shifted 
away from the center of the corresponding pixel and 
the signal in the surrounding pixels must be asym¬ 
metrical. This asymmetry can be removed using the 
expected 10" pointing accuracy if there are no tempo¬ 
ral variations in the Jupiter’s flux or no mirror degra¬ 
dation effect. 

Thus, the map-making algorithm should reflect di¬ 
rectly the scan strategy and the position of Jupiter 
and Saturn (see Fig(T] and |4|. Using the — 30 dB 
threshold, we can estimate the number of pixels, 
which manifests the beam shape in the map for each 
frequency channel. For the simple Gaussian approxi¬ 
mation (see Eq.([2T])) we get 

(22) 

As one can see from Eq. (E^l) . N « 1570 pixels in the 
30 GHz channel (FWHM ~ 33'), while N ss 36 pixels 
only for the 217 GHz (and higher) channels (FWHM 
~ 5'), which is obviously not enough to accurately 
determine these beam shape ellipticity, particularly if 
the ellipticity is no larger than 1.2. 

6. CONCLUSIONS 

In summary, regarding the issues related to using 
transits of planets, such as Jupiter and Saturn, as 
a method of calibration of the in-flight beam shape, 
we conclude the following: 

• The high accuracy of the Cg estimation by 
Planck will require the main beam estimation with 
an error of 1%, which implies that we need to mea¬ 
sure the Jupiter’s flux variation 5B/B = Ss to the 
same level. 

• For observations in the LFI frequency range, 
e.g. 30 ~ 100 GHz, we have a limit of possible varia¬ 
tions of the Jupiter’s flux <0.1 for the —30 dB thresh¬ 
old of the beam. In addition, we need to measure the 
variations of the Jupiter’s flux by using ground-based 
radio telescopes. 
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• When observing at high frequencies 
(2174-857 GHz) we may have problems during obser¬ 
vation of the planets such as missing a target due to 
the narrow beams. 

• Practically, during the mission (ss 15 months), 
we will be able to test the main beam about 3 times 
down to —23.5dB by using the Jupiter’s flux and 
2 times by Saturn. This implies that in the case of 
calibrating the beam degradation effect in intervals 
shorter than 3 months we have to use the method by 
0 . 

• The close (in time interval) transit of the plan¬ 
ets will enable us to check a high-frequency compo¬ 
nent of the beam degradation. 

• Neither Jupiter nor Saturn makes it possible 
to test the far sidelobes. 

• Galactic synchrotron and dust emissions ap¬ 
pear at a level of 0.18 mK. 

• The possible degradation effect could be im¬ 
portant at the same, as the Jupiter’s flux, level of 
variation, ~ 1 ~ 10%, for the beam width. 

• At low frequencies (30 GHz) of the mission, the 
stability of the Jupiter’s flux can be monitored with 
the RATAN-600 radio telescope. Such a possibility 
was demonstrated in investigation of the planet radio 
flux whose stability was checked during a month at a 
level of 0.1% [26] . 

We also note that calibrations by Jupiter and Sat¬ 
urn, together with the method by [7], allow one to 
restore the antenna beam shape for pixelized beam 
on the AT map, which is different from the antenna 
beam shape in the frame on the focal plane. In general 
cases, transition from a pixelized beam to the actual 
beam in the focal plane frame requires a knowledge 
of the noise properties [7J. 
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